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Abstract 



An upper Hmit on the branching ratio for the decay — >■ ti'^vv is set at 
2.4 X 10~^ at the 90% C.L. using pions in the kinematic region 214 MeV/c < 
Pjr < 231 MeV/c. An upper limit of 5.2 x 10^^'' is found on the branching 
ratio for decays ^tt+X", where X° is any massless, weakly interacting 
neutral particle. Limits are also set for cases where Mxo > 0. 
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The process 'w^vV is highly suppressed in the Standard Model by the Glashow- 

Iliopoulos-Maiani (GIM) mechanism [1], which forbids flavor-changing neutral currents 
(FCNC) at tree level. The unequal masses of the virtual up, charm, and top quarks slightly 
spoil the cancellation in one-loop diagrams, allowing a small amplitude to survive. The 
large top quark mass enhances this effect, making the virtual top quark the largest con- 
tributor to ^T^^vV in the Standard Model. Although the charm quark contribution is 
not negligible, recent theoretical calculations [2] have greatly reduced the uncertainty from 
virtual charm quark lines making a measurement of the 'k^vV branching ratio an 

extremely clean measurement of the product of the Cabibbo-Kobayashi-Maskawa quark- 
mixing matrix (Vckm [3]) elements V^gVtd■ Currently, this product is constrained primarily 
by measurements of Vcb and |T46/K;6| through the mixing matrix unitarity [4], but is also 
further restricted by measurements of other higher-order processes, though with significant 
theoretical uncertainty. Combining the recent -^n^uV calculations from Reference [2] 
with constraints on Vckm parameters from CP violation in the neutral kaon system and 
B'^ — B mixing [5], and using the current measurements of the top quark mass [6], the 
K"*" — ^TT+z/I/ branching ratio is most likely in the range (0.6 — 3) x 10^^". New processes such 
as direct FCNC or a decay -^7r~^X^X^, where X° is a new weakly interacting neutral 
particle, could give a three-body branching ratio outside of this range. A two-body decay 
TT'^X'^ would also be new physics, with a cleaner experimental signature. Proposed 
X° candidates are supersymmetric particles [7] or light Goldstone bosons (e.g., Majorons 
[8], familons [9], or axions [10]). 

The E787 detector [11] operating with the low energy separated beam (LESB I) at the 
AGS at Brookhaven National Laboratory collected data for a total of seven months during 
the 1989, 1990 and 1991 running periods. We previously reported 90% confidence level limits 
[12] on the branching ratios, B{K+^Tr+uV) < 5.2 x 10"^ and B{K+^Tr+X°) < 1.7 x 10"^ 
{Mxo = 0) based on data from the 1989 run alone. Here, we add those data to similar 
sets from the 1990 and 1991 runs and analyze the combined set. The 1989-1991 sample 
is the complete data set taken before major upgrades to the beamline and detector were 
undertaken. 

The experiment can be summarized as follows: 800 MeV/c kaons, identified by Cerenkov 
and dE/dX counters, are stopped in a segmented active target (see Figure 1), where they 
decay. Independent measurements of the kinetic energy, momentum, and range (in scintilla- 
tor) of the charged decay products are made using the target, the central drift chamber, and 
the cylindrical range stack of plastic scintillator layers. Pions are distinguished from muons 
kinematically and by identifying the tt —>//—> e decay sequence using 500-MHz transient 
digitizers. Photons are detected in a nearly 47r solid-angle lead-scintillator calorimeter that 
is 12-14 radiation lengths thick. The entire detector is in a 1 T solenoidal magnetic field for 
the momentum measurement. 

The signature for ^tt'^i'V is a decay with a tt"*" as the only observable product. 
The two dominant decay modes are the most important background sources. fJ''^^^ 
{K^2), a two-body decay with a 64% branching ratio, produces a 236 MeV/c K"*" — ^tt+tt" 
{Kt^2), a two-body decay with a 21% branching ratio, produces a 205 MeV/c 7r+. Since the 

^7r~^iyV momentum spectrum extends to 227 MeV/c, we can search for it either above 
or below the K^2 peak. While there is more phase space with < 205 MeV/c, interactions 
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with detector material can shift a ^,^2 pion down into this region making the background 
severe [13]; in this paper, we report on the search above the Kt^2 peak. The dominant source 
of pions with P^^ > 205 MeV/c is beam pions — about two-thirds of the beam particles are 
pions — that scatter from the target into the range stack. 

The search for ^'k^uV follows a threefold strategy: (i) the incident beam particle is 
identified as a that has stopped in the target, with no beam particles at the apparent 
kaon decay time; (ii) the only observed decay product is a single charged-particle track 
identified as a tt"*" that must be delayed in time with respect to the kaon; (iii) the energy, 
range, and momentum of the tt"*" each lie between the Kt^2 and K^2 peaks. A multilevel 
trigger employs each of these elements to reject background events online, while the analysis 
makes more refined use of detector information for further rejection. can survive only 
if the muon is misidentified as a pion and the kinematics are reconstructed incorrectly. 
Kt^2 can survive only if both photons from the 7r° decay are missed and the kinematics are 
reconstructed incorrectly. Scattered beam pions can survive only if the 7r+ is misidentified 
as a with the scattered track mismeasured to be delayed, or if it is missed by the beam 
counters and follows a . The measures taken to deal with the main background sources 
are also very effective against other backgrounds, such as radiative Kf,2 decays or charge 
exchange interactions followed by K^^'w^l^vi. Backgrounds from other decay modes 
were examined and found to be negligible. 

After establishing the overall analysis strategy, we adjusted the cuts with the intention of 
reducing the total expected background to significantly under one event in the final sample. 
The final cuts used were developed during studies of the known background processes. In 
these studies, we take advantage of the redundant methods available for the rejection of each 
background by dividing the cuts used to suppress it into two groups. One group of cuts is 
relaxed or inverted to enhance the background sample, then the other group is applied and 
its rejection is measured. This background-study technique allows us to use data to infer 
background levels of less than one event. For example, a large sample of K^,2 background 
events is obtained by removing the transient digitizer particle identification cuts, and this 
sample is used to measure the K^2 rejection of the kinematic analysis. Similarly, the transient 
digitizer rejection is measured with kinematically selected muons. Assuming these rejections 
are independent, they are combined and used to estimate the number of events that 
will survive the full analysis. Correlations between the two groups of cuts will introduce 
an error in the background estimates from this method, and we group the cuts to minimize 
these effects. 

The detector calibration procedures and analysis software used for the final analysis pre- 
sented here [14] have been refined considerably since Reference [12]. Improvements included 
increased acceptance of the transient digitizer particle identification cuts and improved kine- 
matic resolutions with reduced kinematic tails. An initial analysis [15] had been completed 
before the calibration and software improvements were finished, observing background in 
excess of predictions. The final analysis had a significantly higher acceptance and did not 
suffer from some anomalies in the transient digitizer signals and the kinematic reconstruction 
that may have affected the background predictions in the initial analysis. 

The background from K^,2 (including fi^Vn^y) was evaluated by separately mea- 

suring the rejections of the transient digitizer particle identification and kinematic cuts, and 
is estimated to be less than 0.15 events. The background from Kt^2 was evaluated by sep- 
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arately measuring the rejections of the photon veto and kinematic cuts, and is estimated 
to be less than 0.14 events. The background from beam pion scattering was evahiated by 
separately measuring the rejections of the beam counter and timing cuts, and is estimated 
to be less than 0.07 events. Monte Carlo studies indicated that the background from 
charge exchange interactions was about 0.1 events. 

Figure 2(a) shows the range in scintillator versus kinetic energy for charged tracks in 
the final sample. Only events with a measured charged track momentum in the accepted 
region 211 < < 243 MeV/c are plotted. The rectangular box defines the search region 
in kinetic energy (115 < T^^ < 135 MeV, corresponding to 213 < < 236 MeV/c) and 
range (34 < Ptt < 40 cm of scintillator, corresponding to 214 < P,r < 231 MeV/c), and 
encloses the upper 15% of the 'w^vV phase space. Note that the range cut defines 

the kinematic search region. No events are observed in the signal region. There are six 
events above the signal region, which is consistent with the 5.5 ± 0.6 expected from the Kf,2 
background study. The events clustered at Tj^ = 108 MeV and Pjr = 30 cm are i^^2 decays 
where both photons from the tt'' are missed. The number of such events is consistent with 
Monte Carlo estimates of the photon detection inefficiency [16]. 

Where possible, we used calibration data taken simultaneously with the physics data 
for the acceptance calculation. We relied on Monte Carlo for only the solid angle coverage, 
the accepted region of the 7r+ spectrum, and the losses from tt"*" nuclear interactions and 
decays in fiight. The tt"*" spectrum for ^'k^uV was calculated using a Standard Model 
matrix element with massless neutrinos [17]. Figure 2(b) shows the spectrum of Monte- 
Carlo-simulated T^^vV after the full analysis. Kf,2 calibration data were used to 
measure losses from the beam analysis, from the charged track reconstruction inefficiency, 
from the A'+^7r+ delayed coincidence requirement, and from accidental energy depositions 
at the kaon decay time above our approximately 1 MeV photon veto threshold. Scattered 
beam pion data were used to measure the acceptance of the transient digitizer and kinematic 
TT+Z/i"*" separation cuts. The acceptance calculation is summarized in Table I, resulting in 
a total acceptance of 0.0027 for K+^n+iyV and 0.0127 for K+^n+X^ (Mxo = 0). The 
uncertainty in the acceptance has a negligible effect on limits set with these data. 

During typical running conditions, 3 x 10^ kaons entered the stopping target per 1.5 s 
beam spill. We measure the fraction that decayed at rest in the target to be 0.65 using an 
analysis of Kf,2 data and the well known Kf,2 branching ratio. This normalization to Kf,2 
removes some sources of systematic error from our sensitivity. Our final measured exposure 
for these data is 3.49 x 10^^ stopped kaons. The acceptance (especially the Monte Carlo 
simulation of tt"*" nuclear interactions) and stopping fraction are checked by measuring the 
K^2 branching ratio, for which we obtain 0.205 ± 0.006, where the error includes counting 
statistics and other items not common to the K"*" — ^tt+z/I? analysis. This is consistent with 
the world average, 0.2116 ± 0.0014 [4]. 

Since no events were observed in the signal region (Figure 2(a)), we obtain 90% confidence 
level upper limits oi B{K+ ^ ir+uV) < 2.4 x lO"'^ [18] and B{K+ ^Tr+X°) < 5.2 x 10-^° 
{Mxo = 0) for these data. Figure 3 shows the limit as a function of Mxo for different 
lifetime assumptions, assuming that the photon veto rejects events where the X° decays in 
the spectrometer. 
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FIGURES 
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FIG. 1. Schematic (a) side and (b) end views showing the upper half of the E787 detector. 
C: beam Cerenkov counter; B4: beam hodoscope; I and T: trigger scintillators; RSPC: multiwire 
proportional chambers. 
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FIG. 2. Charged-track range vs. kinetic energy for (a) data and (b) — )■ n^uV Monte 
Carlo for events satisfying the selection criteria (see text) and having measured momentum 
211 < < 243 MeV/c. The rectangular box indicates the search region for 'k'^vV and 

TT+X" {Mxo ~ 0). The horizontal and vertical dashed lines in (b) are the theoretical 
end-points of ^n'^uV in range and energy respectively. 
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FIG. 3. The solid curve gives the 90% confidence level upper limit on the branching ratio for 
-^TT^X^ as a function of M^o. The dashed curves give 90% confidence level upper limits for 
cases where has a finite lifetime. 
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TABLES 



TABLE I. Acceptance factors for K'^^tt'^vv and if+^7r+X° (M^o = 0). Each table entry 
represents the acceptance from a number of related cuts. 



Category 




7r+X° 


Solid angle 


0.43 


0.43 


TT"*" spectrum 


0.15 


0.73 


TT"*" nuclear absorption 


0.53 


0.50 


TT"*" decay in flight 


0.92 


0.92 


— )■ TT"*" delayed coincidence 


0.75 


0.75 


Tr~^ /n'^ kinematics 


0.87 


0.88 


Tr"*" — > /i"*" transient digitizer tagging 


0.41 


0.41 


IJ,^ — )■ e+ transient digitizer tagging 


0.84 


0.84 


Accidental vetoes 


0.67 


0.67 


Beam analysis 


0.84 


0.84 


Reconstruction 


0.69 


0.69 


Net Acceptance 


0.0027 


0.0127 
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